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ABSTRACT 

We have carried out sub- mm ^^C0( J = 3—2) observations of 6 giant molecular 
clouds (GMCs) in the Large Magellanic Cloud (LMC) with the ASTE 10m sub- 
mm telescope at a spatial resolution of 5 pc and very high sensitivity. We have 
identified 32 molecular clumps in the GMCs and revealed significant details of 
the warm and dense molecular gas with n(H2) ~ 10^^^ cm~^ and T^in ~ 60 K. 
These data are combined with ^^C0( J =1 — 0) and ^^C0( J =1 — 0) results and 
compared with LVG calculations. The results indicate that clumps we detected 
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are distributed continuously from cool (~ 10 - 30 K) to warm (~ higher than 30 
- 200 K), and warm clumps are distributed from less dense (~ 10^ cm~^) to dense 
(~ 103-5-5 cm-3).We found that the ratio of ^^CO{J = 3 - 2) to ^^CO{J =1-0) 
emission is sensitive to and is well correlated with the local Ha flux. We interpret 
that differences of clump propeties represent an evolutionary sequence of GMCs in 
terms of density increase leading to star formation. Type I and II GMCs (starless 
GMCs and GMCs with HII regions only, respectively) are at the young phase of 
star formation where density does not yet become high enough to show active star 
formation and Type III GMCs (GMCs with HII regions and young star clusters) 
represents the later phase where the average density is increased and the GMCs 
are forming massive stars. The high kinetic temperature correlated with Ha flux 
suggests that FUV heating is dominant in the molecular gas of the LMC. 

Subject headings: Magellanic Clouds — galaxies: individual (LMC) — ISM: 
clouds — ISM: molecules — radio lines: ISM — submillimeter 



1. Introduction 

It is of a fundamental importance in astronomy to understand the evolution of galaxies. 
Since a major constituent of galaxies is stars, the formation of stars is a fundamental process 
in galactic evolution. The properties of stars characterize the basic contents of galaxies 
and their time evolution. We understand from studies of the Milky Way galaxy that giant 
molecular clouds (GMCs), whose mass ranges from 10^ to 10^ M©, are the principal sites of 
star formation and that it perhaps holds the true in other galaxies as well. We also recognize 
that the CMC properties (e.g., Leo ^ liiie width relation, index of mass spectrum) are similar 



among five galaxies in the Local Group according to the spatially resolved studies (IBlitz et al. 



20061 ). This supports the idea that studies of GMCs will be useful in understanding the 
fundamentals of galactic evolution through the formation and evolution of GMCs and star 
formation therein. 

Observational studies of GMCs have been most effectively made by the mm interstellar 
carbon monoxide emission line at 2.6 mm which allows us to probe molecular gas whose 
density is greater than ~ 100 cm~'^. We note that the most abundant species, molecular 
hydrogen, does not have appropriate line emissions in the mm and sub-mm region due to its 
zero permanent electric dipole moment and large separation between the lowest energy levels, 
which are not excited significantly in the typical physical conditions of molecular clouds. 



Recent advances in sub-mm observations have allowed us to determine physical param- 
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eters of molecular clouds over much larger ranges than in the mm region by comparing line 
intensities between different transitions. These sub-mm studies were initiated by the SEST 
15m telescope in Chile followed by instruments in Hawaii, Mauna Kea and in the Swiss 
Alps at an altitude range from 3700 m to 4200 m, including the CSO 10m, JCMT 15m, 
and KOSMA 3m telescopes, and the AST/RO 1.6m telescope in Antarctica. Subsequently, 
in the 2000 's, the developements of new instruments at an altitude of ~ 5000 m in Ata- 
cama in northern Chile resulted in a superior capability because of the high altitude and dry 
characteristics of the site. The instruments installed in Atacama include the ASTE 10m, 
APEX 12m and NANTEN2 4m telescopes. All these instruments are beginning to take new 
molecular data with significantly better quality than before in terms of noise level as well 
as angular resolution. It is also noteworthy that the current frequency coverage extends as 
high as the 800GHz band and even the THz region. 

Among nearby galaxies we can observe at reasonably high spatial resolutions, the Large 
and Small Magellanic Clouds offer us a unique opportunity to achieve the highest resolutions 
due to their unrivaled closeness, 50 - 60 kpc. In particular, the Large Magellanic Cloud 
(LMC) is actively forming stars in clusters and is an ideal laboratory for us to study star 
formation, particularly massive star formation in star clu ster. In the LMC, t he metallicity 
is a factor of ~ 3 low er than in the Solar neighborhood (jPufour et al.l Il982l : iDufourl Il984j : 



RoUeston et al.ll2002l ). Also, the vis ual extinctions are lower and the FUV field is stronger 
in the LMC than in the Milky Way (jlsrael et al .11 198 61 ). characterizing the initial conditions 
of star formation. 

The first spatially resolved complete sample of GMCs in a single galaxy has been ob- 
tained towards t he whole LMC with the NANTEN 4m telescope in 2.6 mm CO emission at 
40 pc resolution JPukui et al.lll999l . boOll . boO?! : iMizuno et aPboOlh . These studies revealed 
the three types of GMCs in terms of star formation activities; Type I is starless. Type II 
is with HII regions only, and Type HI is associated with active star formation indicated by 
huge HII regions and young star clusters, where the stars identified are only O stars due to 
the sei isitivity limitation. It also revealed t hat t he lifetime of a CMC is as short as ~ 30 
Myrs (IFukuil l2006al : iKawamura et al.l l2006l . 120071 ). These previous studies naturally place 
the LMC as one of the prime targets for sub-mm studies to derive the physical parameters 
of GMCs. 

Another aspect which deserves our attention is that very young, rich stellar clusters are 
forming in the LMC. These are so called populous clusters which are very rare in the Milky 
Way and resemble globulars formed in the primeval Milky Way. The open clusters forming 
in the Milky Way are small in the number of stars and loose in spatial distribution. Along 
with the low metallicity of the LMC, it is an interesting possibility to use molecular data to 
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investigate the formation mechanism of super star clusters at the molecular cloud stage. 



J 



In the past, there have been some studies that used the higher transitions (J = 2 — 1, 
3 — 2, J = 4 — 3, J = 7 — 6) of CO spectra of the mo l ecular clouds in the LMC (e.g.. 



Sorai et al. 


2001 


2003; 


Kim et al. 



2004 



Kim 



Johansson et al.lll998l : iHeikkila et al.lll999l : iBolatto et al.ll2005l : llsrael et al. 

200j. dies suggest the molecular gas may be warmer 



and/or denser than in the Milky Way. 



Johansson et al.l (119981 ) used the SEST 15m telescope to observe the central part of the 
30 Doradus nebula (rms ~ 0.2 K in 0.5 km s~^ velocity resolution for J = 1 — and rms ~ 
1.0 K in 0.5 km s~^ velocity resolution for J = 3 — 2), and the southern HII regions N 158C, 
N 159, and N 160 with a few prominent CO clouds in the J = 2 — 1 and J = 3 — 2 transitions 
of CO. They find that the kinetic temperatures are 10 - 80 K and the highest temperature is 
towards 30 Dor. Th e smallest b earn size a nd grid spacing are 15" and 11" respectively in the 
J = 3 — 2 emission. iHeikkila et al.l (119991 ) used SEST to observe the J = 3 — 2 transition of 
CO in N 159 and 30 Doradus, as well as other rarer molecular species. This study aimed at 
obtaining chemical abundance, while it also provides more information on cloud temperature 
etc, from C0( J = 3 — 2) data. The kinetic temperatures the y derived are 5 K in 30 Dor-10, 
15 K in 30 Dor-27, and 20 - 25 K in N 159W and N 160. ISolatto et all (|2005h employed 
the AST/RO to observe the ^'^CO{J = 4-3) transition at 461 GHz with a 109" beam. 
They observed 9 regions in the LMC at 6'x6' field all with HII regions and derived kinetic 
temperatures from a comparison between the CO(J = 4 — 3) and (J = 1 — 0) transitions. 
N 48, N 55A, N 79, N 83A, N 113, N 159W, N 167, N 214C, and LIRL 648 are included. 
They derive temperatures of 100 - 300 K and note a trend that higher temperatures occur 
in moderate density regions, 100 - 1000 cm^^, and the lower temperatures in much denser 
regions 10^~^ cm^^. These studies were prec eded by a suggest io n that signi f icant amounts of 



warm molecular gas may exist in the LMC (jlsrael et al.ll2003l ). iKim et al.l (120041 ) also made 
similar observ ations towa rds an HII region, N 44, and suggest very dense gas of ~ 10^ cm~^. 
Most recently, iKimI (120061 ) derived Tkin = 100 K, n ~ 10^'^ cm~^ for 30 Dor from the intensity 



ratios of i2C0(J = 7 - 6) to ^^CO{J = 4-3) and ^^CO{J = 1 - 0) to i3C0(J =1-0). 

In the present study, we aim to obtain sub-mm molecular data at better S/N ratios 
than in the previous studies to make estimates of temperatures and densities over a large 
sample in the LMC. We will combine the ^^CO(J = 3 — 2) data obtained with the ASTE 
telescope and CO(J =1 — 0) data obtained with the SEST and Mopra telescopes. In order 
to make reasonable comparisons between the two transitions, J = 3 — 2 and J = 1 — 0, we 
shall smooth the ASTE results (22" beam) to the same resolution as the SEST data (45") 
and use LVG calculations to estimate density and temperature. We shall also employ the 
^^CO(J =1 — 0) data where available to place constraints on the physical parameters. 
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This paper is organized as follows: Section 2 describes the observations. Section 3 and 
4 show the results and data analysis, respectively. In section 5, we discuss the physical 
properties of clumps and evolutional sequence of GMCs. Finally, we present a summary in 
section 6. 



Observations 



2.1. Selection of GMCs 



The pre sent targets were c hosen from the NANTEN catalog of ^^CO(J =1 — 0) GMCs 
compiled by iFukui et al.l (120071 ) . This ca talog is based on th e 2nd survey, with a factor of ~ 
2 higher sensitivity than the 1st survey (IMizuno et al.ll200ll ). The catalog includes 272 CO 
clouds 230 of which are detected at three or more observed positions and they are classified 
into the three Types; 56 (24.3 %) Type I (starless) GMCs, 120 (52.2 %) Type II GMCs (those 
with HII regions only), and 54 (23.5 %) Type III GMCs (those with HII regions and young 
star clusters), where "stars" refer only to O stars due to the limited sensitivity of existiri g 
observations (IKawamura et al.l 120071 : lFukuill2006al : iKawamura et al.ll2006l : iBlitz et al.ll2006l ). 



Among these GMCs, we mainly focus on Type III GMCs whose ^^CO(J = 3 — 2) 
intensities are expected to be strong due to the highly excited conditions. In addition, 
relatively high resolution ^^CO and ^'^CO J = 1 — data, observed with the SEST 15m or 
Mopra 22m telescopes, are collected to derive physical properties of molecular clouds through 
the LVG analysis at a 10 pc scale. 

In the present study, we observe GMCs in the south-east region of the LMC, which 
contains 30 Doradus - the largest and most massive HII region in the Local Group. We 
observe the molecular ridg e extending southw ard from 30 Doradus, and " CO Arc" along the 
southeastern optical edge JFukui et al.lll999h . 3 Type III GMCs, LMC N J0538-6904 (the 
30 Dor region), LMC N J0540-7008 (the N 159 region and the N 171 region), and LMC N 
J0530-7106 (the N 206 region), are selected as the principal targets, and two Type II GMCs, 
LMC N J0544-6923 (the N 166 region) and LMC N J0532-7114 (the N 206D region), and 
a Type I CMC, LMC N J0547-7041 (the CMC 225 region), are included for reference. The 
locations of the observed GMCs and regions are shown in Figure [H and their coordinates 
and the data used in this paper are summarized in Table [H Hereafter, the region names, 
which are in the parenthesis above or column 4 of Table [H are used to identify the regions. 
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2.2. i2cO(J = 3-2) 



Observations of the ^^CO(J = 3 — 2) transition at 345GHz were made with th e ASTE 
(Atac ama Submilhmeter Telescope Experiment) telescope at Pampa la Bola in Chile (lEzawa et al 
2004 ) in October 2004. The half power beam width was measured to be 22" at 345 GHz 
by observing the planets. This corresponds to 5.3 pc at the distance of the LMC, 50 kpc. 
The telescope was equipped with a single "cartridge-type" SIS receiver , sensitive from 324 
to 384 GHz, which is of a similar design t o that for ALMA (jKohndl2005l ). The spectrometer 
was an XF-type digital auto-correlator (ISorai et al.l l2000l ). and was used in the wideband 
mode, which has a bandwidth of 512 MHz with 1024 channels. The spectrometer provided 
a velocity coverage and resolution of 450 km s~^ and 0.44 km s~^, respectively, at 345 GHz. 
We observed 6 GMCs (7 regions) in the Large Magellanic Cloud as shown in Figured] and 
listed in Table [TJ These observations were carried out by position switching at a grid spacing 
of 20" or 30" for the entire clouds, and of 10" or 15" for the regions around the local peaks of 
the integrated intensity. The pointing error was measured to be within 7" in peak to peak by 
observing CO point sources R Dor or o Get every 2 hours during this observing term. The 
spectral intensities were calibrated by employing the standard room-temperature chopper- 
wheel technique. We observed Ori-KL once a day, and N 159W every 2 hours to check the 
stability of intensity calibration, and the intensity variation during these observations was 
less than 13 %. We use 0.7 for the main-beam efficiency at 345 GHz, which was measured 
by observing Jupiter. The system noise temperature was typically 300 K in double-sideband 
(DSB) including the atmosphere towards the zenith. The typical r.m.s. noise fluctuations 
were ~ 0.25 K at a velocity resolution of 0.44 km s~^ for a 1 minute integration for an 
on-position. In total, about 1400 points were observed in Equatorial coordinates (B1950). 
Velocities were relative to the Local Standard of Rest (LSR). These observations were made 
remotely from an ASTE operation room of San Pedro d e Atacama, Chile, usin g the network 
observation system N-C0SM0S3 developed by NAOJ (IKamazaki et al.ll2005l ). 



2.3. i2co(J = 1-0) and i3C0(J= 1-0) 

2.3.1. Mopra observations 

A 20' X 120' region, the prominent molecular ridge extending from 30 Doradus south- 
ward, was mapped in the J = 1 — transition of ^^CO at a frequency of 115 GHz with the 
22m ATNF Mopra telescope, in 5 runs from May to September 2005. This region contains 
the 30 Dor, N 159 and N 171 regions. The newly implemented on-the-fiy (OTF) mode was 
used, in which the telescope takes data continuously while moving across the sky. Spectra 
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were taken at a 6" spacing so that the 33" Mopra beam would be well oversampled in the 
scanning direction; the row spacing was 8", also assuring oversampling. The typical system 
noise temperature, Tgys, was 500 K in the single side band (SSB) towards the zenith. The 
pointing was checked on the SiO maser R Dor every 2 hours; typical pointing error was less 
than 5" rms. The digital correlator was configured to output 1024 channels across 64 MHz 
in each of two orthogonal polarizations. The velocity resolution and coverage were 0.16 and 
160 km s~^, respectively at 115 GHz. The observing time was about 100 minutes per field 
(5'x5'), providing rms noise fiuctuations of ~ 0.34 K at a velocity resolution of 0.65 km s^^. 
Initial spectral processing (baseline fitting and calibration onto a scale) was performed 
using the livedata task in AIPS++, and the spectra were gridded into data cubes using the 
AIPS++ gridzilla task. During the gridding, a Gaussian smoothing kernel with a FWHM 
taken at a half of the beam size was convolved with the data, so the effective resolution of 
the output cubes was 36" . The cubes were then rescaled onto a scale using an "extended 
beam" efficiency of 0.55 (ILadd et al.ll2005l ). This takes into account that sources larger than 
about 80" in diameter will couple to both the main beam and the inner error beam of the 
telescope. The cubes were 2 dimensional Gaussian smoothed to a 45" beam, which is the 
beam size of SEST at 115 GHz. 



2.3.2. SEST observations 

Observations towards N 206, N 206D and GMC 225 regions in the ^'^CO{J =1-0) line 
(115 GHz) were made in August, 2001 and February, 2002, using the SEST 15m telescope 
at La Silla, Chile. The HPBW was 45" at 115 GHz, the front end was the IRAM 115 SIS 
receiver and the spectrometer was a high-resolution AOS with 2048 channels. The typical 
system noise temperature was 550 K (SSB). The velocity resolution and coverage were 0.2 
and 216 km s~^ respectively at 115 GHz. We mapped these 3 regions in position switching 
with a grid spacing of 40" or 20". The typical integration time was 1 minute for an on- 
position, providing rms noise fiuctuations of ~ 0.18 K at a velocity resolution of 0.2 km 

Observations towards N 206, N 206D and GMC 225 regions in the ^^CO{J =1-0) line 
(110 GHz) were made in February and December, 2002, also using the SEST 15m telescope. 
We mapped peak positions of ^^CO(J = 1 — 0) in position switching with a grid spacing of 
20". The system noise temperature was typically 230 K (SSB). The typical integration time 
was 4 minutes for an on-position, providing rms noise fiuctuations of ~ 0.04 K at a velocity 
resolution of 0.2 km s~^. 

The pointing accuracy was 5" rms. We checked this by observing SiO maser toward R 
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Dor every 2 hours during this observing term. N 159W was observed periodically for pointing 
checks and intensity calibration. We use 0.8 for the main-beam efficiency at 115 GHz by 
assuming main-beam tem perature T^h of N 159W to be ~ 6.5 - 6.9 K to keep consistency 
with former publications ([Johansson et al.lll994l . Il998l ). 



Observations toward N 166 region are described separately by iGaray et al.l (120021 ). 



3. Results 

We first present the ^^CO(J = 3 — 2) images of the clouds at 5 pc resolution and make 
an empirical comparison of them with the ^^CO(J =1 — 0) distribution (section 3.1.). Next, 
we define molecular clumps and estimate the physical parameters of each clump (section 
3.2.). 



3.1. Distributions of the ^^CO{J = 3-2) emission 

In Figure ^ typical ^^CO{J = 3-2) and ^^CO{J = 1-0) profiles of the 30 Dor, 
N 159, and GMC 225 regions are presented. The upper panels. Figures [2]^a) - (c), show 
^^CO(J = 3 — 2) profiles. These illustrate the low noise levels of the present data, typically, 
~ 0.25 K rms at 0.44 km s~^ velocity resolution. Among the present observed positions, 
the ^^CO{J = 3-2) intensity is strongest at T^i, ~ 12.3 K towards N 159W (Figure [2](b) 
). The lower panels. Figures [21(d) - (f), show ^^CO(J = 1 — 0) profiles towards the same 
positions, where the red line indicates ^^CO(J = 3 — 2) profiles convolved to a 45" Gaussian 
beam following the method described in Section 4. The ^^CO(J = 3 — 2) intensities from 
the 30 Dor and N 159 regions are a little stronger than the ^^CO(J = 1 — 0) intensity when 
convolved to the same resolution. Only towards GMC 225 is the ^^CO(J = 3 — 2) intensity 
about 50% weaker than the ^^C0( J = 1—0) intensity. The peak velocity and line width are 
nearly the same between the two transitions in these regions. 

The distributions of the integrated intensities of ^^C0( J = 3 — 2) and ^^C0( J = 1 — 0) 
are shown in Figures [31 -[HI Detailed descriptions of each region are presented in the following. 



3.1.1. 30 Dor (Figure\^ 



Figures[la) and (b) show the distributions of the ^'^CO{J = 3-2) and ^'^CO{J =1-0) 
integrated intensities in the 30 Dor region. We see a general trend that the J = 3 — 2 
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distribution shows more details which are not obvious in the J = 1 — distribution owing to 
the higher angular resolution and possibly due to the more compact distribution of warmer 
and denser gas in J = 3 — 2 than in J = 1 — 0. We see three peaks corr esponding to the 



12 



C0( J = 1 — 0) peaks, 30 Dor-10, 30 Dor-6, and 30 Dor-12, reported in [Johansson et al. 



(119981 ). One of them in the north, 30 Dor-06, which is singly peaked in J = 1 — appears 
to be resolved into two peaks with the present beam. 



3.1.2. N 159 (Figure\^ 

FiguresSKa) and (b) show the ^^CO{J = 3-2) and ^'^CO{J =1-0) distributions in N 
159. We note that N 159W shows the strongest intensity among the present clouds as well 
as a very compact peak which is not sufficiently resolved with the present beam. Its radius 
is estimated to be a few pc after de-convolution. N 159E also shows a compact distribution 
with a hint of a sub peak, while N 159S shows similar distributions both in J = 1 — and 
J = 3 — 2. The east-west elongation of N 159S may be due to the scanning effect of the OTF 
mapping, and needs to be confirmed. 



3.1.3. N 171 (Figure\§ 



Figures[5](a) and (b) show the ^^GO{J = 3-2) and ^^CO{J =1-0) distributions in N 
171. We note that the J = 3 — 2 emission is weaker than the J = 1 — emission. There are 
multiple velocity components, at Vlsr = 225 km s~^, 230 km s ~^, and 240 k r n s~^ in both 
J = 3 — 2 and J = 1 — as indicated in Table [2] and Table 1 of iKutner et al.l (119971 ) . 



3.1.4. N 166 (Figure\^ 



Figures El^a) and (b) show the ^^CO{J = 3-2) and ^^GO{J =1-0) distributions in 
N 166. The J = 1 — data themselves have already been published by Garay et al. (2002). 
We see four peaks in both J = 3 — 2 and J = 1 — 0. Of these peaks, three peaks are named 
Cloud-B, Cloud-C, and Cloud-D, respectively, as reported by I Garay et al.l (120021 ). although 
Cloud-C is resolved into two peaks with the present beam and observing grid. 



3.1.5. N 206 (Figure^ 



Figures El^a) and (b) show the ^'^CO{J = 3-2) and ^'^CO{J =1-0) distributions in 
N 206. A J = 1-0 peak appears to be resolved into two sub peaks and a north-south 
filamentary structure with the higher angular resolution of the J = 3 — 2 line. 

3.1.6. N 206D (FigureW^ 

FigureslHKa) and (b) show the ^^00(7 = 3-2) and ^^00(7 =1-0) distributions in N 
206D. We can see a head-tail structure in J = 3 — 2, although it appears more rounded in 
J = 1 - 0. 

3.1.7. GMC 225 (Figure \B) 

Figures Ha) and (b) show the ^'^CO{J = 3-2) and ^'^CO{J =1-0) distributions in 
GMC 225. We note that J = 3 — 2 emissions are weaker than J = 1 — 0. 

3.2. Properties of the clumps 

We identified clumps according to the following way in the J = 3 — 2 distributions shown 
in Figures [3] - [HI 1) Pick up local peaks using the integrated intensity. 2) Draw a contour at 
half of the peak integrated intensity level and identify it as a clump unless it contains other 
local peaks. 3) When there are other local peaks inside the contour, draw new contours at 
the 70 % level of each integrated intensity peak. Then, identify clumps separately if their 
contours do not contain another local peak (the boundary is taken at the "valley" between 
clumps), or else identify a clump by using the highest contour as a clump boundary. 4) 
If a spectrum has multiple velocity components with a separation of more than 6 km s~^, 
identify these components to be associated with different clumps. 

As a result, 32 clumps have been identified; 5 clumps in 30 Dor region, 10 clumps in N 
159 region, 6 clumps in N 171 region, 5 clumps in N 166 region, 2 clumps in N 206 region, 
1 clump in N 206D region, and 3 clumps in GMC 225 region. Their line parameters at the 
peak positions are shown in Table [2l and their physical properties are listed in Table O The 
clump size, line width, and virial mass range 1.1 - 12.4 pc, 4.0 - 12.8 km s~^, and 4.6x10^ - 
2.2x10^ Mq, respectively. The smallest clump identified by the procedure above is detected 
with 2 observing points and its size is 1.1 pc. There are several other weak emissions (below 
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6 a), and they are detected with 1 observing point. These emissions are not identified as 
clumps. When we assume such small clumps (i? ~ 1 pc, dV ~ 7 km s~^) are exist, their 
virial mass is 9.3 x 10^ Mq, and this corresponds to detection limit. There are 18 local peaks 
above the 6 a noise level which are not identified as clumps, because the stronger peaks are 
located near these local peaks. Their virial mass is 6.5 x 10^ M0, if their size and line width 
are similar to identified clumps (i? ~ 7 pc, dV ~ 7 km s~^). This seems to correspond to 
completeness limit. 

If the clumps are identified with the 70 % level of the peak integrated intensity level 
instead of the 50 %, the clump size changes to a half of the original one, althogh the line 
width does not change. Their virial mass also changes to about a half of the original one. 
Histograms of their physical properties are presented in Figure [101 Typical values are 7 pc, 

7 km s-\ and 6x10^ M©, in size, line width, and virial mass, respectively. Thei r line widths 



are larger than those of the GMCs in our galaxy (e.g., IWilliams fc Blitall998l : llkeda et al. 



I999I : ISun et al.ll2006l ). We note a trend that the masses of these clumps are relatively large 
compare d with those o f high mass cloud cores in th e Milky Way which are of the order of 10'^ 



M0(e.g., iBurton et al.ll2005l : lYonekura et al.ll2005l ). Precise comparisons will be necessary 



using same tracers and resolution in the future, because those galactic studies above are 
based on optically thin mm dust continuum and C^^O emission. The clumps in the LMC 
are also fairly compact, with sizes of several pc or less. Hereafter, the region names and the 
numbers of clumps are used to identify clumps (e.g., "30 Dor No.l"). 



4. Data Analysis 

4.1. Derivation of line intensity ratios 

The spatial resolution of the present CO data varies depending on the telescope and 
frequency. We convolved and regridded these data into the same resolution and position 
using 2 dimensional Gaussian smoothing in order to compare them to derive reliable peak 
intensity ratios of ^^CO{J = 3 - 2) to ^^CO{J = 1-0) (hereafter, i?3_2/i_o). 

We made Gaussian fits to each of the ^^CO(J = 3 — 2) and ^^CO(J =1—0) spectra 
having a single peak in most cases. We derived peak intensities and FWHM line widths 
through these fittings. The ratios of the two transitions were then derived as the ratio 
between peak TmbS. The distributions of these ratios are shown along with the ^^CO(J = 
1 — 0) distributions in Figures [TT] - [T71 (a) . The youngest stellar clust ers SWBO, whos e ages 



are estimated to be less than 10 Myrs, are also shown by red circles (IBica et al.l 119961 ). 



Averaged i?3_2/i-o over each clump (hereafter, R3-2/i-o,ciump) was also derived from the 
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averaged ^^C0( J = 3 — 2) and ^^CO(J = 1 — 0) spectra over the each clump. A summary 
of the R3-2/i-o,dump for each clump is shown in Table HI The histogram in Figure [18] shows 
that the R3-2/i-o,dump ranges from 0.2 to 1.6. These ratios will be compared with numerical 
calculations of radiative transfer in the LVG approximation to derive constraints on density 
and temperature. 



4.2. LVG analysis 

4-2.1. Calculations of LVG model 

To estim ate the physical propertie s of the molecular gas in the LMC, we performed an 
LVG analysis ( Goldreich fc KwanI 1974 ) of the CO rotational transitions. The LVG radiative 
transfer code simulates a spherically symmetric cloud of constant density and temperature 
with a spherically symmetric velocity di stribution pro portional to the radius, and employs 
a Castor's escape probability formalism (jCastoii 119701 ). It solves the equations of statistical 
equilibrium for the fractional population of CO rotational levels at each density and tem- 
perature. It includes the lowest 40 rotational levels of the ground vibra tional level and uses 
the Einstein's A and H2 impact rate coefficients of ISchoier et al.l (120051 ). 



The present calculations incorporate the lowest 40 rotational levels of CO in the ground 
vibrational state over a kinetic temperature range of T^in = 5 - 200 K and a density range 
of n{B.2) = 10 - 10^ cm~^. We did not include higher energy levels in the present study, 
which requires including populations in the lower vibrationally excited states. Therefore, 
the present work does not cover kinetic temperatures above 200 K, which should be dealt 
with in the future for analyses of higher sub-mm transitions above J = 4 — 3. This imposes 
a limit of Tkin ~ 200 K in the present study and even higher temperature is not excluded in 
general below. 

We performed calculations for 3 different CO fractional abundances; V(CO) = [C0]/[Il2] 
= lxlO-^ 3x10"^ and lxlO-^ and 3 different ^^CO/^^cQ abundance ratios of 20, 25, and 
30 (Heikkila et al. 1999). 



4.2.2. Results from ^^CO(J = 3 - 2) and ^^CO(J = 1 - O; data for 32 clumps 

First, we assume that V(CO) is uniform among the clumps and derive density and 
temperature using the LVG results. 



Figure [TW a) illustrates that we obtain the following lower limits for kinetic temperature 
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and density. For the 8 clumps with R3-2/i-o,ciump > 1 (30 Dor No.l, 2, 3, 4, N 159 No.l, 2, 
6, 8), we estimate T^in >30 K and n(B.2) >10'^ cm~^. For the 24 clumps with R3-2/i-o,dump 
<1 (30 Dor No.5, N 159 No.3, 4, 5, 7, 9, 10, N 171 No.l, 2, 3, 4, 5, 6, N 166 No.l, 2, 3, 4, 5, 
N 206 No.l, 2, N 206D No.l, GMC 225 No.l, 2, 3), we estimate T^n >several K and ^(Ha) 
>severalxlO^ cm~^. In either case, the lower limits tend to increase with the R3-2/i-o,ciump- 

4.2.3. Results from CO(J = 3 - 2;/^ C0(J = l-Q), and CO(J = 1 - O; data for 13 

clumps 

We can better constrain these physical parameters when ^^CO(J = 1 — 0) data are 
available. Figure [19] shows the general behavior of the loci of constant /23_2/i-o and constant 
i?i2/i3 (peak intensity ratio of ^^C0( J = 1—0) to ^^C0( J = 1—0)) in the density-temperature 
plane. It is recognized that the combination of the two will allow us to determine the 
parameters relatively well since the two lines are nearly "orthogonal" in the plane, except 
for densities higher than ~ 10*^ cm~^ (Figure [T^ . 

Of the 32 clumps, ^^C0( J =1 — 0) data are available for 13 clumps, including 4 clumps 
of R3-2/ 1-0, dump > 1 (30 Dor No.l, 4, N 159 No.l, 2) and 9 clumps with R3-2/i-o,dump <1 
(N 159 N0.4, N 166 No.l, 3, 4, N 206 No.l, 2, N 206D No.l, GMC 225 No.l, 3). For these 
13 clumps, we made a detailed analysis using the -R12/13 at peaks of ^^CO(J = 1 — 0) and 
have obtained the best constraints. 

We summarize the input parameters for the 13 clumps in columns 3-5 of Table [51 All 
the data refer to the ^^CO(J =1 — 0) beam size, 45"; the higher-transition data has been 
Gaussian smoothed as described in section 4.1. Clump averaged dv/dr were used for the 
calculations. R3-2/i-o,dump and -R12/13 at the peak of ^^CO(J =1—0) were used, and the 
errors of these ratios are both estimated as ±20 %, which are derived from errors of absolute 
intensity calibration. These correspond 27 a and 7 a noise levels of ^^CO(J = 3 — 2) and 
^^CO(J = 1 — 0), respectively. This indicates the errors of intensity ratios are dominated by 
the error of absolute intensity calibration. The way of the clump definition does not change 
R3-2/ 1-0, dump, but changes dv/dr, and dv/dr is about 2 times as large as the original one, 
when the clumps are identified with the 70 % level of the peak integrated intensity level. 
This, however, does not affect the results of the LVG calculations. 

We describe three typical cases; (a) 30 Dor No.l, (b) N 206 No.l, and (c) GMC 225 
No.l are shown in Figure [22] A fractional CO abundance of X{CO) of 3xl0~^ was used 
throughout. The horizontal axis is molecular hydrogen density (n(H2)), and the vertical axis 
is the gas kinetic temperature (Tkm). Solid lines are R3-2/i-o,dump and dashed lines are i?i2/i3 
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at ^^CO(J = 1 — 0) peaks. Hatched areas indicate the overlap regions of these two ratios 
within the errors which is allowed from the observed ratios. It also includes the uncertainty 
due to a possible variation of the ^^CO/^'^CO abundance ratio from 20 to 30. Hereafter, we 
shall call clumps above 30 K "warm" and those below "cold". Clumps with densities greater 
than 10^'^ cm~^ are referred to as "dense" and those at lower densities as "less dense". 

30 Dor No.l exhibits a R3^2/i-o,dump of ~ 1.4, suggesting that it is a warm (Tkin= 50 
- 200 K) and dense {n = lO^"^ cm^^) clump (Figure [2n](a)). N 206 No. 1 (Figure [2D](b)) and 
CMC 225 No.l (Figure I^UT c)) show of R3-2/i-o,ciump ~ 0.7 and ~0.4, respectively. Figures 
[20](b) and (c) indicate that they are warm (Tkm = 30 - 200 K) and less dense clumps (n = 
10^ cm~^) and cold (Tkm = 15-40 K) and less dense (n = 10^ cm~^) clumps, respectively. 
We find that -R12/13 is useful to discriminate the temperature difference except in the case of 
30 Dor No.l. 

All the figures for the 13 clumps for all different fractional CO abundances (1x10^^, 
3xl0~^, and 1x10^^) are given in Appendix A. In Table [5] we present the results of for the 
13 clumps for a fixed fractional CO abundance of X(CO) = 3x10^^. The results of previous 
studies are also summarized in Table for comparison. We find that the present results 
are the most extensive among these studies in terms of the number of samples, while they 
are basically consistent with previous study for the individual clumps or clouds. There are 
relatively large differences in both density and temperat ure in N 159 No.4 (N 159S) only. 



and this is because of the high density tracers are used in iHeikkila et al.l (119991 ) . 



Figure [2T] summarizes the calculated densities and temperatures for the 13 clumps. We 
see first that the temperature ranges from 10 K up to more than 200 K and density from 10'^ 
cm~^ to 10^ cm~^. Clumps we detected are distributed continuously from cool (~ 10 - 30 
K) to warm (~ higher than 30 - 200 K), and warm clumps are distributed from less dense 
(~ 10^ cm~^) to dense (~ lO^-^~^ cm~^), although cool clumps are all less dense. The three 
cases shown in Figure [20] represent typical cases, "warm and dense" (Tkm ~ higher than 30 
- 200 K, n{B.2) ~ lO^-^-^ cm'^), "warm and less dense" (Tkin ~ higher than 30 - 200 K, 
n(H2) - 10^ cm-3), and "cool and less dense" (Tkm ~ 10 - 30 K, ^(Ha) ~ 10^ cm'^). 



4.2.4. Effects of X (CO) 

We will now discuss the possible effect of changing X(CO), as summarized in Table 
[61 If we adopt X(CO) = lxlO~^, the contours of -R3-2/1-0 shift to lower temperature 
and the contours of -R12/13 shift to higher density. Accordingly, the solution shifts to lower 
temperature, higher density. If we adopt X{GO) = lxlO~^, the i?3_2/i-o contours shift to 
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higher temperature and the -R12/13 contours shift to lower density. Accordingly, the solution 
shifts to higher temperature, lower density (Table [6]). 



Next, we vary X{CO) from clump to clump. iHeikkila et al.l (119991 ) estimate X{CO) in 
30 Dor No.l (30 Dor-10), N 159 No.l (N 159W), and N 159 No.4 (N 159S) to be 1x10"^ 
lxlO~^, and 3xl0~^, respectively. Table [6] indicates that 30 Dor No.l shows similar values 



Heikkila et al. lJl999[). wh ereas N 159 No.l and No.4 become warmer and lower in density 



to 



than in 



Heikkila et al. 



19991). This d iscrepancy may be due to the fact that the high density 



tracers used by iHeikkila et al.l ( 1l999l ) are not used in the present study. To summarize, the 
assumption of uniform fractional abundance, X{CO), is fairly good and the present results 
do not show significant difference even if we adopt the different fractional abundances used 
by previous authors. 



4.3. Comparisons to Ha flux 



4-3.1. Relation between R3-2/1-0 o,nd Ha flux 



We converted the Ha data (IKim et al.l Il999l ) towards the present clouds using the 



10 ergs 



-10 gj^gg g 1 (^j^ 



^ towards strong HIT regions 
3 — 2) data grids. Ha flux images 



method given in the Appendix B. The typical background level of the Ha flux is 
s~^ cm~^ at the 40" scale and ranges up to 10 
(Figure [22I). These data were re-gridded i nto the ^^COf J 

with the youngest stellar clusters SWBO (IBica et al.lll996l . younger than 10 Myrs) are shown 
in Figures fTTl - [T71 (b) and the ^^CO(J = 3 — 2) contours from Figures [3] - [9] are overlayed 
for comparison in these figures. 

Figures [TT] - [T71 indicate a clear trend that the i?3_2/i-o is enhanced to 1.0 - 1.5 towards 
HH regions or clouds with young clusters, as in 30 Dor and N 159 regions. On the other 
hand, the ratio is low, less than 1.0, towards clumps with neither HH regions nor clusters, 
as in GMC 225. We also note that the ratio is enhanced towards the regions where Ha 
is intense or towards the interfaces between clouds and HH regions. A summary of the 
averaged Ha flux over each clumps is shown in Table |H It is not fully guaranteed that 
the Ha emission is actually in contact with the molecular gas and some of the apparent 
coincidence could be fortuitous. Nonetheless, previous studies comparing CO and Ha flux 
indicate a strong correlation between them and lend a su pport to the assumption that almost 



all the coincidences indicate actual physical association (IFukui et al.lll999 
200l[ iKawamura et al.lboOTl l 



Yamaguchi et al. 



Figure 1231 shows the correlation between R-i-2/i~o,ciump and the averaged Ha flux over 
each clump. It is clear that i?3-2/i-o,dnmp is well correlated with the averaged Ha flux. 
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with a correlation coefficient of ~ 0.79. This is a fairly good empirical relationship and 
should be tested to see if it holds true in other galaxies. The relation suggests that higher 
i?3__2/i-o reflect higher temperatures or higher densities. It is notable that the clumps with 
the averaged Ha flux greater than 10~^^ ergs s~^ cm~^ correspond to high i?3-2/i-o,d«mp 
of 1.0 - 1.5. In the barred spiral galaxy M83, the CO(J = 3 - 2)/(J = 1 - 0) integrated 
intensity ratio exceeds unity at the nucleus, whereas the ratio gradually de creases to 0.6 



O/Twith distance from center. The ratio is constant through the disk region (IMuraoka et al. 
20o3). 



4-3.2. Comparison with physical properties and Ha emisson 

Figure [2l] shows plots of density and temperature as functions of averaged Ha flux, and 
gives us another insight into these properties. Where averaged Ha is strong (10~^^ergs s^^ 
cm~^) the clumps are always warm at around T^in = 100 K or more. On the other hand, 
when averaged Ha is weak, density is always low but temperature can be either high or low. 

The high density clumps show high i?3_2/i-o,dump of 1.0 - 1.5 and are associated with 
strong Ha flux, while the low density clumps show low R3-2/i-o,dump of 0.5 - 1.0 with weak 
Ha flux. 

Since we are averaging the Ha intensity in each clump, we may be diluting the localized 
low Ha flux towards some of the clumps in Type HI and II GMCs. The effects of Ha emission 
or nearby clusters on the molecular gas are perhaps local phenomena as indicated by the 
comparison in Figures [TT] - [T7I for the individual clumps. 



5. Discussions 

5.1. Dense and compact clumps as candidates for proto-cluster condensations 

We have carried out sub-mm ^^C0( J = 3 — 2) observations of GMCs in the LMC which 
are most extensive and highly sensitive compared to the previous studies. Six GMCs were 
selected based on the NANTEN CO survey of the LMC, including 3 Type HI GMCs actively 
forming O stars in addition to 3 Type I/II GMCs which are quiet in 0-star formation or 
cluster formation, although the formation of low to intermediate mass star is not excluded. 
The spatial resolution of ~ 5 pc and the high sensitivity allowed us to identify 32 molecular 
clumps in these GMCs and to reveal significant details of the warm and dense molecular gas 
with n(H2) ~ 10^-5 cm'^ and Thh ~ 10 - 200 K. 
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The typical mass of the molecular clumps is large, in the range of 5x10^ - 2x10^ 
Mq with radii of 1 - 12 pc. Of all of our objects, N 159 No.l or -W shows the strongest 
concentration of mass of ~ 7x10^ Mq within a radius of ~ 5 pc. The masses seem to be 
larger than those of typical Milky Way GMCs such as those in the eta Car region (e.g.. 



Yonekura et al.ll2005l ). although the propeties of these galactic GMCs are based on optically 



thin C^^O data. We suggest that these are good candidates for the precursors of rich super 
clusters like R136 in 30 Dor which includes more than 10^ st ars in a small volume with a 
radius of ~ 1 pc. It is of particular interest to look for even denser gas towards them in 
higher excitation transitions of the sub-mm region. 



5.2. Density and temperature of the clumps and implications 

The results of our LVG analysis indicate that clumps are distributed from cool to warm 
in temperature and from less dense to dense in density. These differences of clump properties 
in density and temperature show good correspondence with the CMC Types based on the 
star formation activity, as well as with the Ha emission of ionized gas associated with each 
clump. Clumps in Type III GMCs are warm (Tkm ~ 30 - 200 K) and are either dense (n(H2) 
~ 10^'^"^ cm~'^) or less dense {n{B.2) ~ lO^cm"'^). Clumps in Type II GMCs are either warm 
(Tkin ~ 30 - 200 K) or cool (Tkm ~ 10 - 30 K) and less dense (^(Hs) ~ 10^ cm'^). Clumps 
in Type I CMC are cool (Tkin ~ 10 - 30 K) and less dense {n{B.2) ~ lO'^cm"^). The physical 
parameters of clumps are generally correlated with the star formation activity of GMCs and 
can perhaps be interpreted in terms of evolutionary effects. 

Our interpretation is that defferences of clump density and temperature represent an 
evolutionary sequence of GMCs in terms of density increase leading to star formation; Type 
I/II GMCs are at a young phase of star formation where density has not yet reached high 
enough values to cause active massive star formation, and Type III GMCs represent the later 
phase where the average density is higher, including both high and low density sub-types. 
The high density clumps in Type III GMCs show high R3-2/i-o,dump of 1.0 - 1.5 and are 
associated with strong Ha flux while the low density clumps in Type HI GMCs show low 
Rz~2/i-o,ciump of 0.5 - 1.0 with weak Ha flux. 

We suggest two alternative ideas to explain the density difference of the clumps in Type 
HI GMCs; one is that density is being enhanced by shock compression driven by HII regions 
and the other is that gravitational condensation of each clump plays a role in the density 
increase. The former may be difficult because the shock front may occupy a small volume 
which is likely missed with the present 5 pc beam. It seems thus favorable that the latter 
scenario is working mainly to enhance density. 
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The present study, which resolved the smaller clumps in GMCs at 5 - 10 pc scales, 
indicates that the clumps may have physical properties affected by local properties such 
as the Ha distribution. It should be interesting to investigate the variations among these 
internal clumps and their relation to star formation. 



5.3. FUV Heating of the molecular gas in the LMC 



The present findings that the R3-2/i~o,dump is well correlated with Ha flux suggests that 
the heating of molecular gas by far-ultraviolet (FUV) photons may be effective in the LMC 
where the dust opacity is lower and the FUV intensity is higher than in the Milky Way. 
The molecular gas in the Milky Way is mainly heated by cosmic ray protons of ~ 100 MeV 
as discussed by a number of authors, although the surface layers of molecular clouds with 
small visual extinctions at Av ~ a few mag or less may be dominated by the FUV heating 



[e.g., iKaufman et al.lll999l ). Some authors have ma de detailed calculation s of gas heating 
and cooling under the effects of FUV radiation fileds ( Kaufman et al.l 1999 ). We shall try to 
present a picture that can be applied to the present results below. 

First, the gas temperature is d etermined through the balan ce between the cooling and 
heating. According to Table 4 in [Goldsmith fc Langerl (119781 ) . the total cooling rate is 
6.8x10-27^2-2 ergs cm'^ s"i for X (CO) / {dv / dr) = 4x10"^ and ^(Hs) = 10^ cm'^. In 30 
Dor region, since X(CO)/( dv/dr) = 3xl0~^ / 0.8 = 3 .75 xlO~^ and this value is 10 times 
lower than the value used in [Goldsmith fc Langen (119781 ). n(H2) = 10'^ cm~'^ can be read 10'^ 



cm ^, then the cooling rate is estimated as 1.7x10 22 ergs cm ^ s ^ for T 
value is a factor 2-3 smaller than that of Galactic clouds with n{ll2) = 10^ 



cm 



100 K. This 

3 and T = 



50 K. 



We shall assume that the heating by cosmic ray electrons is not important in the LMC. 
This assumption is not directly confirmed, but it is consistent with the low non-thermal 
fraction of the LMC's radio continuun emission (e.g., iHughes et al.ll2006l ) and studies that 



suggest a sign ificant fra c tion of cosmic ray elect rons are able to escape from low luminosity 
galaxies (e.g.. [Bei[2003l : Iskillman fc Kleinlll988h 



The FUV flux 



1995 



Go) is estimated as 3500 for 30 Doradus ( Bolatto et al. 



Werner et a. 



Israel et al. 



19961: 



1978 



1999 



Israel fc Koornneej ll979h. and 300 for N 159 JSolatto et al 



Poglitsch et al 



Israel fc Koornneef 19791 ). In Orion, it is estimated as 25 ( Bolatto et al. 



1999 



I999I : Istacev et allll993[ ). The FUV flux in the LMC is larger than that in the Milky Way. 



PDR models are calculated by iKaufman et al.l (119991 ) which incorporate the chemical 
and physical processes that form and destruct atoms or molecules, as well as ionization 
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effects. Figuer 1 of iKaufman et al.l (Il999l ) shows the kinetic temperature for a molecular 
gas layer with density of n (cm~^) under FUV flux of Gq at the surface. PDR surface 
temperatures are estimated as listed in Table [71 These indicate that temperature becomes 
as high as 100 - 300 K on the PDR surface under the conditions of the clumps in Type III 
GMCs in the LMC These temperatures are basically consistent with the temperatures of 
the warm clumps in the present sample. 

Generally speaking, at a scale of ~ 10 pc, Tkin ~ 100 K seems to be higher than the 
kinetic temperatures typical in Milky Way GMCs, where the Milky Way v alues are usually 



derived from the ^^CO(J = 1 — 0) emission only (e.g., eta Car Tkm ~ 50 K by lYonekura et al. 



20051 ). This suggests that the heating of molecular clouds may be stronger in the LMC than 
in the Milky Way and the molecular temperature may be higher. If this is correct, the lower 
metallicity, resulting in lower extinction, is the basic cause for the higher temperature in 
addition to the stronger FUV field in the LMC. We shall note in the end that this higher 
temperature in the molecular gas possibly leads to an increase of the Jeans mass of molecular 
clumps, which may favor the formation of rich super clusters in the LMC. This is consistent 
with the higher mass of the molecular clumps which may represent precursors of the clusters. 

The present work has undertaken to sample 6 GMCs (7 regions) to have a uniform 
determination of the density and temperature in the LMC. The number of GMCs is still 
limited to 6 among ~ 300 detected with NANTEN. We should make more efforts to collect 
appropriate data sets in the sub-mm wavelengths to improve our understanding of the cloud 
properties. NANTEN2, ASTE and others will certainly be powerful in achieving this goal. 



6. Summary 

We summarize the results as follows. 

1) We have used the ASTE 10m telescope to obtain the distribution of ^^C0( J = 3 — 2) 
emission at 345 GHz towards 6 GMCs (7 regions) in the LMC. We have identified 32 clumps 
in these GMCs at ~ 5pc resolution. The radius, line width and virial mass are estimated as 
1.1 - 12.4 pc (7 pc), 4.0 - 12.8 km s"^ (7 km s^^), and 4.6x10^ - 2.2x10^ Mq(6x10^ Mq), 
respectively, with the average values in the parenthesis. 

2) We have compared the present results with LVG radiative line transfer calculations in 
order to obtain the density and temperature estimated for the 13 clumps using R3^2/i-o,ciump 
and -R12/13 at 45" resolution. The clumps are distributed from cool (~ 10 - 30 K) to warm 
(more than ~ 30 - 200 K) and from less dense (n(H2) ~ 10^ cm~^) to dense (n(Il2) ~ 
cm~^) 
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3) The Ha flux towards these clumps is well correlated with the ^^C0( J = 3—2) /^^C0( J = 
1 — 0) ratio, R3-2/i-o,dump, and clumps with Ha fluxes greater than 10~^^ ergs s~^ cm~^ have 
large R3-2/i-o,dump of ~ The ^^C0( J = 1 — 0) data were taken with the SEST 15m and 
Mopra 22m telescopes. 

4) The typical mass of the molecular clumps ranges 5x10^ - 2x10^ M© with radii of 1 
- 12 pc. Of all of our objects, N 159 No.l or -W shows the strongest concentration of mass 
of ~ 7x10'^ M0 within a radius of ~ 5 pc. We suggest that these are good candidates for 
the precursors of rich super clusters like R136 in 30 Dor. 

5) We suggest that differences of clump properties represent an evolutionary seqeunce 
of GMCs in terms of density increase leading to star formations. Type l/ll GMCs are at 
a young phase of star formation where density has not yet reached high enough values to 
cause active massive star formation, and Type HI GMCs represent the later phase where the 
average density is higher, including both high and low density sub-types. 

6) The high kinetic temperature correlated with Ha flux suggests that FUV heating is 
dominant in the molecular gas of the LMC. The low fraction of non-thermal radio continuum 
emission and calculations of PDR models support this suggestion. Furthermore, the high 
temperature in the molecular gas possibly leads to an increase of the Jeans mass of molecular 
clumps, which may favor the formation of rich super clusters. 
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A. Complete LVG results of all clumps 

Figures [231 - 1571 show the complete LVG results described in section 4.2. We present all 
cases: X{CO) of (a)lxlO-^ (b)3xl0-^ and (c)lxl0-5 for 13 clumps. 
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B. Ha flux 



We shall he re describe the method of scaling data values in FITS cube (Ha image in 
Kim et al.l (119991 ): hereafter, "Kim's FITS"). The data values in Kim's FITS is not flux scale, 
and the calibration is needed for quantitative comparison with CO clouds and i?3_2/i-o- The 
procedu re is as follows. 1) Sum up the data values of Kim's FITS inside apertures which are 
listed i nlKennicutt fc Hodgd (119861) for each hsted HII region. 2) Plot cataloged values of Ha 
flux in iKennicutt fc Hodgd (119861 ) as a function of summed values derived in 1). They are 
well fitted by a power function oi y = 5.16 x 10^^^x°'^ (c.c. = 0.94) (Figure 1381) . 3) Convert 
data values of Kim's FITS to Ha flux scale (ergs s"^ cm^^) using the function derived above. 



C. LVG results in the other planes 

C.l. Physical properties - i?3-2/i-o,d«mp (Figure [39]) 

The density plots (Figure [5W a) and (b)) show that higher -R3-2/i-o,d«mp (>l-0) corre- 
spond to higher densities of 10^ to 10^ cm~^, while lower Rz-2/i-o,ciump (<l-0) correspond 
to lower densities of around 10^ cm~'^. The temperature plots (Figure ISW c) and (d)) show 
that higher i?3-2/i-o,c/«mp (>0.5) correspond to higher temperatures of >30K, while lower 
Rz-2/i-o, clump (<0.5) correspond to lower temperatures of <30K. 

Then, roughly speaking, we can say that clumps with i?3-2/i-o,dnmp lower than 0.5 
have lower densities of around 10^ cm~^ and lower temperatures of <30K, clumps with 
R3^2/ 1-0, dump of 0.5 to 1.0 havc lower density around 10^ cm~^ and higher temperatures of 
>30K, and clumps with Rz-2/i-Q,ciump higher than 1.0 have higher densities of 10"^ to 10^ 
cm~^ and higher temperatures of >30K, although ratios, densities, and temperatures are 
with large error bars. 

Of course, there are great benefits to using i?i2/i3 in LVG analyses. We could not obtain 
the above results with only i?3_2/i_o, as mentioned in section 4. 



C.2. Physical properties — i?i2/i3 (Figure 1401) 

The density (Figure HOlfa) and (b)) does not show a significant correlation with i?i2/i3. 
The temperature (Figure 1401(c) and (d)) shows a good correlation with i?i2/i3, that is, higher 
ratios indicate higher temperatures. Usually, i?i2/i3 correspond to density, but in this case, 
due to the LVG analysis using both i?3_2/i_o and -R12/13, larger -R12/13 indicate lower density 
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and lower density tends to higher temperature. 



C.3. i?i2/i3 - Ha flux (Figure E]) 

There is not significant relation between i?i2/i3 and Ha flux, contrary to -R3_2/i-o ratio. 
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Fig. 1. — CO velocity integrated intensity map (iFukui et al.l l! 
image. The observed 7 regions are indicated by white circle. 
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Fig. 2.— i2Q0(J = 3-2) and ^^00(7 =1-0) profiles at selected points. Note that 
the vertical scale is main beam temperature, Tmb. The velocities (abscissae) are relative 
to the LSR. (a)i2C0(J = 3-2) profile of 30 Dor. {hy^CO{J = 3-2) profile of N 159. 
(c)i2C0(J = 3-2) profile of GMC 225. (d) ^^CO{J =1-0) and smoothed ^^CO{J = 3-2) 
profiles of 30 Dor. (e) ^^CO{J = 1-0) and smoothed ^^CO{J = 3-2) profiles of N 159. 
(f) i2C0(J =1-0) and smoothed ^^CO{J = 3-2) profiles of GMC 225. 
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Fig. 3. — (a) Contour map of ^^CO(J = 3 — 2) integrated intensity in the 30 Doradus region. 
The contour levels are 5, 10, 15, 20, 25, 30, 35, 40, and 45 K km s~^ Observed points 
are indicated by dots, (b) Contour map of ^^CO(J = 1—0) integrated intensity in the 30 
Doradus region. The contour levels are 10, 15, 20, and 25 K km s~^. 
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Fig. 4. — (a) Contour map of ^^C0( J = 3 — 2) integrated intensity in the N 159 region. The 
contour levels are 5, 10, 20, 30, 40, 50, 60, 70, 80, 90 and 100 K km s'^ Observed points 
are indicated by dots, (b) Contour map of ^^CO(J = 1 — 0) integrated intensity in the N 
159 region. The contour levels are 10, 20, 30, and 40 K km s~^. 
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Fig. 5. — (a) Contour map of ^^CO(J = 3 — 2) integrated intensity in the N 171 region. 
The contour levels are 8, 12, 16, and 20 K km s~^. Observed points are indicated by dots, 
(b) Contour map of ^^CO(J =1—0) integrated intensity in the N 171 region. The contour 
levels are 8, 12, 16, 20, 24, 28, 32, 36, 40, 44, 48, and 52 K km s-\ 
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Fig. 6. — (a) Contour map of ^^CO(J = 3 — 2) integrated intensity in the N 166 region. 
The contour levels are 6, 8, 10, 12, 14, 16, 18, 20, and 22 K km s~^. Observed points are 
indicated by dots, (b) Contour map of ^^CO(J = 1 — 0) integrated intensity in the N 166 
region. The contour levels are 8, 10, 12, 14, 16, 18, 20, and 22 K km s~^. Observed points 
are indicated by dots. 
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Rigtit Afecanssion (B1 950) Righrt MwnMiqri ^;B1 950) 

Fig. 7. — (a) Contour map of ^^C0( J = 3 — 2) integrated intensity in the N 206 region. The 
contour levels are 4, 6, 8, 10, 12, 14, 16, 18, 20, 22, 24, and 26 K km s'^ Observed points 
are indicated by dots, (b) Contour map of ^^CO(J = 1 — 0) integrated intensity in the N 
206 region. The contour levels are 4, 6, 8, 12, 14, 16, and 18 K km s^^. Observed points are 
indicated by dots. 



RighlAitsnsian (B195(J) RIglit Aw»nslon (B1S&0) 



Fig. 8. — (a) Contour map of i2C0(J = 3-2) integrated intensity in the N 206D region. 
The contour levels are 4, 6, 8, 10, 12, 14, 16, and 18 K km s~^. Observed points are indicated 
by dots, (b) Contour map of ^^CO(J =1 — 0) integrated intensity in the N 206D region. 
The contour levels are 4, 6, 8, 10, 12, 14, 16, 18, 20, 22, 24, and 26 K km s'^ Observed 
points are indicated by dots. 
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Fig. 9. — (a) Contour map of ^^CO(J = 3 - 2) integrated intensity in the CMC 225 region. 
The contour levels are 3, 4, 5, 6, 7, and 8 K km s~^. Observed points are indicated by dots, 
(b) Contour map of ^^CO(J =1 — 0) integrated intensity in the CMC 225 region. The 
contour levels are 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, and 19 K km s~\ 
Observed points are indicated by dots. 
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Fig. 10. — Histograms of physical properties of ^^CO(J = 3 — 2) clumps: (a)size, (b)line 
width, and (c)virial mass. 
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Fig. 11. — (a) Color map of -R3_2/i-0; and (b) Ha flux image of the 30 Doradus region. 
Contours are ^^CO(J = 3 — 2) integrated intensity. Contour levels are the same as Figure 
3a. Thick lines indicate the observed area in ^^CO(J = 3 — 2), and the red circle indicates 
the position of young cluster (<10 Myrs; SWBO). 
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Fig. 12. — (a) Color map of -R3_2/i-o, and (b) Ha flux image of the N 159 region. Contours 
are ^^CO(J = 3 — 2) integrated intensity. Contour levels are the same as Figure 4a. Thick 
lines indicate observed area of ^^CO(J = 3 — 2), and the red circles indicate the position of 
young clusters (<10 Myrs; SWBO). 
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Fig. 13. — (a) Color map of -R3_2/i_o, and (b) Ha flux image of the N 171 region. Contours 
are ^^CO(J = 3 — 2) integrated intensity. Contour levels are the same as Figure 5a. Thick 
lines indicate observed of i2C0(J = 3-2). 
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Fig. 14. — (a) Color map of -R3_2/i-o, and (b) Ha flux image of the N 166 region. Contours 
are ^^CO(J = 3 — 2) integrated intensity. Contour levels are the same as Figure 6a. Thick 
lines indicate observed area of ^^CO(J = 3 — 2). 
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Fig. 15. — (a) Color map of i?3_2/i_o, and (b) Ha flux image of the N 206 region. Contours 
are ^^CO(J = 3 — 2) integrated intensity. Contour levels are the same as Figure 7a. Thick 
lines indicate observed area of ^^CO(J = 3 — 2). 
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Fig. 16. — (a) Color map of i?3_2/i-o, and (b) Ha flux image of the N 206D region. Contours 
are ^^CO(J = 3 — 2) integrated intensity. Contour levels are the same as Figure 8a. Thick 
lines indicate observed area of ^^CO(J = 3 — 2). 



Fig. 17. — (a) Color map of i?3_2/i-o, and (b) Ha flux image of the CMC 225 region. 
Contours are ^^CO(J = 3 — 2) integrated intensity. Contour levels are the same as Figure 
9a. Thick lines indicate observed area of ^^CO(J = 3 — 2). 
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Fig. 18. — Histgram of clump averaged peak intensity ratio of ^^C0( J = 3 — 2) to ^^C0( J 

1 — 0) {R3-2/l-0,dump) ■ 
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Fig. 19. — Contour plots of LVG analysis for reference. Contours are (a) -R3-2/1-O7 (b) -R12/13, 
and (c) a+b. X(CO)=3xlO-^ dt;/dr=1.0km s'^ pc-\ and abundance ratio of i^CO/^^co 
is 25. 
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Fig. 20.— Contour plots of LVG analysis of 3 clumps: (a)30 Dor No.l, (b)N 206 No.l, and 
(c)GMC 225 No.l. The vertical axis is kinetic temperature T]^in, and the horizontal axis is 
molecular hydrogen density n{B.2)- Solid lines are R3-2/i-o,dump, and dashed lines are -Ri2/i3- 
Hatched areas are the regions in which these two ratios overlap within intensity calibration 
errors of 20% and uncertainty due to a possible variation of ^^CO/^'^CO abandance ratio 
from 20 to 30. 
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Fig. 21. — Plot of LVG results. The vertical axis is kinetic temperature, Tkm, and the 
horizontal axis is molecular hydrogen density n{B.2). 
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Fig. 22.- 
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Histgram of clump averaged Ha flux. Background level is nearly 10 ergs s ^ 
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Fig. 23.— Plots of i?3 -2/1-0, dump as a function of clump averaged Ha flux by (a)region and 
(b)GMC type, respectively. 
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g. 24. — Plots of physical properties as a function of clump averaged Ha flux, (a) 
region, (b) ^(Ha) by GMC type, (c) Tkin by region, (d) Tkm by GMC type. 
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Fig. 25. — Contour plots of LVG analysis of 30 Dor No.l. The vertical axis is kinetic 
temperature, Tkm, and the horizontal axis is molecular hydrogen density, n{B.2). Solid lines 
are R3-2/i-o,dump, and dashed hues are Ru/is- Hatched areas are the regions in which 
these two ratios oberlap within intensity calibration errors and uncertainty due to a possible 
variation of ^^qo/i^cO abandance ratio. X{CO) = (a)lxlO-^ (b)3xl0-^ and (c)lxlO-^ 
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Fig. 26. — Contour plots of LVG analysis of 30 Dor No. 4. The vertical axis is kinetic 
temperature, Tkm, and the horizontal axis is molecular hydrogen density, n(H2). Solid lines 
are R3^2/i-o,ciump, and dashed hues are Ru/is- Hatched areas are the regions in which 
these two ratios oberlap within intensity calibration errors and uncertainty due to a possible 
variation of ^^cQ/iscO abandance ratio. X(CO) = (a)lxlO-^ (b)3xl0-^ and (c)lxlO-^ 
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Fig. 27. — Contour plots of LVG analysis of N 159 No.l. The vertical axis is kinetic tem- 
perature, Tkin, and the horizontal axis is molecular hydrogen density, n{B.2). Solid lines are 
R3-2/ 1-0, clump, and dashed lines are -Ri2/i3- Hatched areas are the regions in which these two 
ratios oberlap within intensity calibration errors and uncertainty due to a possible variation 
of i^CO/i^CO abandance ratio. X(CO) = (a)lxlO-^ (b)3xl0-^ and (c)lxlO-^ 
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Fig. 28. — Contour plots of LVG analysis of N 159 No. 2. The vertical axis is kinetic tem- 
perature, Tkin, and the horizontal axis is molecular hydrogen density, n(H2). Solid lines are 
R3-2/ 1-0, dump, and dashed lines are -Ri2/i3- Hatched areas are the regions in which these two 
ratios oberlap within intensity calibration errors and uncertainty due to a possible variation 
of i^CO/i^CO abandance ratio. X(CO) = (a)lxl0-6, (b)3xl0-^ and (c)lxlO-^ 
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Fig. 29. — Contour plots of LVG analysis of N 159 No. 4. The vertical axis is kinetic tem- 
perature, Tkin, and the horizontal axis is molecular hydrogen density, n(H2). Solid lines are 
Rz-2/i-o,ciumpi and dashed lines are -Ri2/i3- Hatched areas are the regions in which these two 
ratios oberlap within intensity calibration errors and uncertainty due to a possible variation 
of i^CO/i^CO abandance ratio. X(CO) = (a)lxl0-6, (b)3xl0-^ and (c)lxlO-^ 
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Fig. 30. — Contour plots of LVG analysis of N 166 No.l. The vertical axis is kinetic tem- 
perature, Tkin, and the horizontal axis is molecular hydrogen density, n(H2). Solid lines are 
i?3_2/i-o,dMmp, and dashed lines are -Ri2/i3- Hatched areas are the regions in which these two 
ratios oberlap within intensity calibration errors and uncertainty due to a possible variation 
of i^CO/i^CO abandance ratio. X(CO) = (a)lxl0-6, (b)3xl0-^ and (c)lxlO-^ 
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Fig. 31. — Contour plots of LVG analysis of N 166 No. 3. The vertical axis is kinetic tem- 
perature, Tkin, and the horizontal axis is molecular hydrogen density, n(H2). Solid lines are 
R3-2/ 1-0, clump, and dashed lines are i?i2/i3- Hatched areas are the regions in which these two 
ratios oberlap within intensity calibration errors and uncertainty due to a possible variation 
of abandance ratio. X(CO) = (a)lxlO-^ (b)3xl0-^ and (c)lxlO-^ 
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Fig. 32. — Contour plots of LVG analysis of N 166 No. 4. The vertical axis is kinetic tem- 
perature, Tkin, and the horizontal axis is molecular hydrogen density, n(H2). Solid lines are 
i?3_2/i-o,cZump5 and dashed lines are i?i2/i3- Hatched areas are the regions in which these two 
ratios oberlap within intensity calibration errors and uncertainty due to a possible variation 
of abandance ratio. X(CO) = (a)lxlO-^ (b)3xl0-^ and (c)lxlO-^ 
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Fig. 33. — Contour plots of LVG analysis of N 206 No.l. The vertical axis is kinetic tem- 
perature, Tkin, and the horizontal axis is molecular hydrogen density, n{B.2). Solid lines are 
R3-2/ 1-0, clump, and dashed lines are -Ri2/i3- Hatched areas are the regions in which these two 
ratios oberlap within intensity calibration errors and uncertainty due to a possible variation 
of i^CO/i^CO abandance ratio. X(CO) = (a)lxlO-^ (b)3xl0-^ and (c)lxlO-^ 
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Fig. 34. — Contour plots of LVG analysis of N 206 No. 2. The vertical axis is kinetic tem- 
perature, Tkin, and the horizontal axis is molecular hydrogen density, n(H2). Solid lines are 
R3-2/ 1-0, clump, and dashed lines are i?i2/i3- Hatched areas are the regions in which these two 
ratios oberlap within intensity calibration errors and uncertainty due to a possible variation 
of abandance ratio. X(CO) = (a)lxlO-^ (b)3xl0-^ and (c)lxlO-^ 
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Fig. 35. — Contour plots of LVG analysis of N 206D No.l. The vertical axis is kinetic 
temperature, Tkin, and the horizontal axis is molecular hydrogen density, n(B.2)- Solid lines 
are R3-2/i-o,dump, and dashed hues are Ru/is- Hatched areas are the regions in which 
these two ratios oberlap within intensity calibration errors and uncertainty due to a possible 
variation of ^^cQ/iscO abandance ratio. X(CO) = (a)lxlO-^ (b)3xl0-^ and (c)lxlO-^ 
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Fig. 36. — Contour plots of LVG analysis of GMC 225 No.l. The vertical axis is kinetic 
temperature, Tkm, and the horizontal axis is molecular hydrogen density, n{B.2). Solid lines 
are -R3-2/i-o,dump, and dashed hues are R^/is- Hatched areas are the regions in which 
these two ratios oberlap within intensity calibration errors and uncertainty due to a possible 
variation of ^^cQ/iscO abandance ratio. X(CO) = (a)lxlO-^ (b)3xl0-^ and (c)lxlO-^ 
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Fig. 37. — Contour plots of LVG analysis of GMC 225 No. 3. The vertical axis is kinetic 
temperature, T^in, and the horizontal axis is molecular hydrogen density, n{B.2). Solid lines 
are R3-2/i-o,dump, and dashed hues are Ru/is- Hatched areas are the regions in which 
these two ratios oberlap within intensity calibration errors and uncertainty due to a possible 
variation of ^^cQ/iscO abandance ratio. X{CO) = (a)lxlO-^ (b)3xl0-^ and (c)lxlO-^ 
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Fig. 38. — Plot of relation between integrated Ha flux and FITS value. The solid line 
represents the least-squares fit. 
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. 39. — Plots of physical properties as a function of R3-2 /i-o, dump- 
n(H2) by GMC type, (c) Tki„by region, (d) T^nby GMC type. 



(a) n(H2) by region. 
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Fig. 40. — Plots of physical properties as a function of -Ri2/i3- (a) n(B.2) by region, (b) 
?T.(H2) by GMC type, (c) Tkmby region, (d) Tkinby GMC type. 
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Fig. 41. — Plots of i?i2/i3 as a function of clump averaged Ha flux (a) by region and (b) by 
GMC type. 



Table 1. List of observed GMCs and transitions 





GMC 








Position 






Telescope 




Num. ^ 


Name ^ 


Type^ 


Region Name 


a(1950) 


(5(1950) 


Ref." 


^^CO(J=3-2) 


^^CO(J=1-0) 


^^CO(J=1-0) 










(h m s) 


(O , „) 










186 


LMC N J0538-6904 


III 


30 Dor 


5 39 0.2 


-69 8 0.0 


3 


ASTE 


MOPRA 


SEST 3 


197 


LMC N J0540-7008 


III 


N 159 


5 40 18.2 


-69 47 0.0 


3 


ASTE 


MOPRA 


SEST 3 








N 171 


5 40 24.1 


-70 8 0.0 


4 


ASTE 


MOPRA 




216 


LMC N J0544-6923 


II 


N 166 


5 44 52.5 


-69 26 39.1 


5 


ASTE 


SEST 5 


SEST 5 


153 


LMC N J0530-7106 


III 


N 206 


5 31 33.9 


-71 10 0.0 


6 


ASTE 


SEST 


SEST 


156 


LMC N J0532-7114 


II 


N 206D 


5 32 52.2 


-71 16 0.0 


6 


ASTE 


SEST 


SEST 


225 


LMC N J0547-7041 


I 


GMC 225 


5 48 35.7 


-70 40 0.0 


6 


ASTE 


SEST 


SEST 



Note. — Column (1); Running number of GMC used in Table 1 in Pukui et al.(2006b). Column (2): Name of GMC. Column (3): Type of GMC. 
Column (4): Region name used in this paper. Column (5)-(7): Coordinates used as reference position in each region for these ^^CO(J = 3 — 2) 
observations. Column (8)-(10): Telescope used for each observation. 

^These coordinates are used as reference positions in each regions for this ^•^CO(J = 3 — 2) observations 

''References for the positions. 

References. — (1) Pukui et al. 2006b; (2) Kawamura et al. 2006; (3) Johansson et al. 1998; (4) Kutner et al. 1997; (5) Garay et al. 2002; (6) 
Mizuno et al. 2001 
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Table 2. Observed properties of ^^CO(J = 3 — 2) clumps 



Position Peak properties 



Region 


No. 


«(1950) 


<5(1950) 




Vlsr 


Ay 


I.I. 


Another Ident. 






(h m s) 


i° 


' ") 


(K) 


(km s ) 


(km s ) 


(K km s^-') 




30 Dor 


I 


5 39 8 fi 


-69 


6 15 


5.2 


250.4 


8.4 


49.5 


30nor-10 1 




2 


5 38 54 6 


-69 


8 


4.3 


246.9 


6.9 


31.7 


30Dor-12 ^ 




3 


5 38 49 


-69 


4 30 


3.1 


251.9 


5.4 


22.5 






4 


5 38 49.0 


-69 


3 30 


2.6 


248.2 


5.0 


21.1 


30Dor-06 ^ 




5 


5 38 54.6 


-69 


9 


1.2 


247.4 


4.5 


5.3 




N 159 


1 


5 40 3.7 


-69 


47 


12.3 


237.3 


8.1 


106.6 


N159-W 1 




2 


5 40 35 5 


-69 


46 


7.6 


233.2 


7.0 


55.5 


N159-E ^ 




3 


5 40 47.1 


-69 


46 15 


5.4 


234.6 


6.7 


46.0 






4 


5 40 32.7 


-69 


52 


3.5 


234.4 


9.7 


35.4 


N159-S 1 




5 


5 39 49.3 


-69 


47 


3.0 


233.5 


9.5 


29.7 






6 


5 40 47.2 


-69 


51 30 


3.6 


235.6 


5.3 


19.8 






7 


5 39 55.1 


-69 


45 


4.5 


238.3 


3.4 


16.6 






8 


5 40 29.8 


-69 


50 


0.91 


234.6 


6.9 


10.5 






9 


5 40 41.4 


-69 


51 


0.87 


235.7 


8.0 


9.3 






10 


5 39 37.7 


-69 


45 30 


0.89 


238.6 


7.9 


8.3 




N 171 


1 


5 40 0.5 


-70 12 30 


2.5 


224.3 


5.2 


15.4 


30DOR-CENTER-04 ^ 




2 


5 40 0.5 


-70 10 30 


1.3 


227.3 


8.3 


12.7 






3 


5 40 24.1 


-70 


9 45 


1.3 


230.9 


9.2 


18.2 






4 


5 40 24.1 


-70 


8 30 


2.3 


231.0 


6.7 


16.1 


30DOR-CENTER-03 ^ 




5 


5 40 6.4 


-70 


9 


1.7 


240.6 


4.9 


11.5 






6 


5 40 18.2 


-70 


9 30 


1.9 


240.7 


5.1 


10.4 




N 166 


1 


5 44 58.2 


-69 26 39 


4.5 


227.6 


4.5 


22.1 


Cloud-C 3 




2 


5 44 46.8 


-69 27 24 


2.9 


229.4 


6.2 


20.3 






3 


5 45 3.9 


-69 29 9 


2.2 


227.9 


5.5 


14.6 


Cloud-D ^ 




4 


5 44 32.6 


-69 


23 39 


1.0 


236.1 


9.2 


10.1 


Cloud-B 3 




5 


5 45 9.6 


-69 


30 9 


0.82 


231.6 


4.5 


9.6 




N 206 


1 


5 31 29.8 


-71 


9 40 


5.1 


229.0 


4.6 


26.2 






2 


5 31 46.3 


-71 


8 20 


2.3 


231.2 


4.6 


13.5 




N 206D 


1 


5 32 58.4 


-71 15 20 


4.5 


224.6 


4.0 


18.9 




GMC 225 


1 


5 47 51.3 


-70 41 20 


1.8 


216.5 


3.9 


9.0 






2 


5 48 25.6 


-70 41 20 


1.2 


216.2 


6.2 


8.6 






3 


5 48 53.9 


-70 40 30 


1.8 


216.4 


3.8 


7.9 





Note. — Column (1): Region. Column (2): Running number in eadi region. Column (3)-(4): Positions of observed points 
with maximum ^^C0( J = 3 — 2) integrated intensities within the ^■^CO( J = 3 — 2) clumps arc given in equatorial coordinates. 
Column (5)-(8): Observed properties of the ^■^CO(J = 3 — 2) spectra obtained at the peak positions of the ^■^CO(J = 3 — 2) 
clumps. The peak main-beam temperature Tj^b, l^SR) ^iid the FWHM line width AV are derived from a single Gaussian 
fitting, and arc given in columns 5, 6, and 7, respectively. In column 8, the ^^CO(,7 = 3 — 2) integrated intensities at the 
peak positions of the ^^CO(J = 3 — 2) clumps are shown. Column (9): Another indentifications based on ^"^CO(J = 1 — 0) 
observations with BEST. 

References. — (1) Johansson et al. 1998; (2) Kutner et al. 1997; (3) Garay et al. 2002 
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Table 3. Physical properties of ^^C0( J = 3 — 2) clumps 



Clump properties 



Region 


No. 


chimp 


-^doconv 








(Icin s 


(1>'-) 




30 Dor 


1 


12.8 


7.2 


22 




2 


6.9 


3.1 


2.9 




3 


5.8 


3.1 


2.0 




4 


4.0 


3.7 


1.1 




5 




1.1 


0.46 


N 159 


1 


8.7 


4.7 


6.8 




2 


7.3 


6.9 


6.9 




3 


7.6 


4.7 


5.2 




4 


8.5 


11.3 


15 




5 


10.2 


4.7 


n A 
9.4 




6 


5.2 


4.3 


2.2 




/ 


A 1 
4.1 


O Q 
Z.O 


U. / 




Q 
O 


D.O 


D.O 


A 1 




n 


V 1 









lU 


' " ' 


2.3 


Z. ( 


N 171 


1 


6.5 


6.6 


5.2 




2 


6.9 


8.3 


7.4 




3 


7.2 


7.8 


7.7 




4 


6.2 


5.9 


4.3 




5 


6.1 


7.5 


5.3 




6 


5.6 


5.9 


3.5 


N 166 


1 


6.1 


9.9 


6.9 




2 


6.5 


9.5 


7.6 




3 


6.3 


10.3 


7.7 




4 


9.2 


12.4 


20 




5 


9.5 


4.7 


8.1 


N 206 


1 


6.6 


7.1 


5.9 




2 


5.2 


6.3 


3.2 


N 206D 


1 


4.0 


7.3 


2.2 


GMC 225 


1 


4.2 


8.8 


3.1 




2 


6.6 


5.5 


4.6 




3 


4.9 


8.8 


3.9 



Note. — Column (1): Region. Column (2): Running num- 
ber in each region. Column (3): Line width, AV^iu^p, de- 
rived by using a single Gaussian fitting for a spectrum ob- 
tained by averaging all the spectra within a single clump. 
Column (4): Deconvolved effective radii, -Rdcconvi defined as 
[-^nodcconv ^ HPBW /2)^] "'^ . .Rnodcconv is effective radii de- 
fined as (A/tt)*^'^, where A is the total cloud surface area. 
Column (5): Virial masses of the i2cO(J = 3-2) clumps, a 
density profile of p oc r~^ was assumed; the virial mass is writ- 
ten as Mvir / M0 = 190 X [AKiumpCkm s^i)]^ x i?deconv(pc) 
(MacLaren et al. 1998) 
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12 


CC)( T — 




CO(J 


= 1 


— 0) ratio and Ha flux 


Region 


INO. 


Tmb(3-2) 


rj-i /-I (W 


^3-2/1-0, clump 


Averaged Hq flux 






(K) 


(K) 








(xlO""*^*^ ergs s~"^ cm""^) 


30 Dor 


1 


1.63 


1.20 




1.4 




56.49 




2 


1.76 


1.12 




1.6 




73.37 




3 


1.56 


1.13 




1.2 




24.09 




4 


1.98 


1.58 




1.3 




12.83 




5 


0.73 


0.87 




0.8 




39.71 


N 159 


1 


5.50 


5.04 




1.1 




9.28 




2 


3.69 


3.34 




1.1 




9.37 




3 


2.92 


3.25 




0.9 




6.36 




4 


2.33 


3.34 




0.7 




1.14 




5 


1.84 


2.19 




0.8 




1.93 




6 


1.61 


1.57 




l.U 




1.42 




7 


1.67 


1.98 




0.9 




3.14 




8 


0.78 


0.57 




1 A 

1.4 




1.52 




9 


0.80 


0.99 




0.8 




1.26 




10 


0.57 


0.84 




0.7 




1.42 


N 171 


1 


1.04 


2.20 




0.5 




1.00 




2 


0.69 


1.55 




0.4 




0.93 




3 


1.02 


3.09 




0.3 




0.97 




4 


1.34 


3.89 




0.3 




0.98 




5 


0.88 


1.92 




0.4 




0.95 




6 


0.87 


2.72 




0.3 




0.97 


N 166 


1 


1.97 


3.09 




0.6 




1.05 




2 


1.59 


2.60 




0.6 




1.01 




3 


1.17 


1.92 




0.6 




1.06 




-1 


0.79 


1.26 




0.6 




1..31 




5 


0.1)2 


0.92 




0.7 




1.01) 


N 206 


1 


1.83 


2.39 




0.8 




3.26 




2 


1.34 


2.44 




0.5 




2.68 


N 206D 


1 


2.26 


4.70 




0.5 




1.11 


GMC 225 


1 


1.04 


2.81 




0.4 




0.93 




2 


0.55 


1.62 




0.3 




0.91 




3 


0.77 


2.09 




0.4 




0.92 



Note. — Column (1): Region. Column (2): Running number in ea<;h region. Column (3); The 
pecik main-beam temperature, T^ji,, of the ^^CO(J = 3 — 2) spectra derived by using a single 
Gaussian fitting for a spectrum obtained by averaging all the spectra to the beam within a single 
clump. The intensities refer to the ^^CO{J = 1 — 0) beam size (45"). Column (4): The peak 
main-beam temperature, Tju^i of the ^^CO( J = 1 — 0) spectra derived by using a single Gaussian 
fitting for a spectrum obtained by averaging all the spectra within a single clump. Column (5): 
The ratios of Tnib(3-2) to Tnib(l-O). Column (6): The Ha flux obtained by averaging within a 
single clump. 



Table 5. Summary of LVG and MPE analyses 



Region 


No. 


dv/dr 


^3-2/1-0, ciump 




L2/13 


n(ti2j (cm 




(P: 


revious 


studies) 






(km s pc ) 




Value 


Ref. 


This worli 


(1) 


(4) 


(5) 


(6) (7) (8) 


30 Dor 


1 


0.9 


1.4 


11.5 


(I) 


3x10* - 3x10® 
>50 


1Q4.5 

40 - 80 


1x10® 
50 




lO** 
100 




4 


0.5 


1.3 


17.7 


(1) 


1x10* - 1x10® 
>60 


10*® 
>20 








N 159 


1 


0.9 


1.1 


8.6 


(1) 


3x10* - 8x10® 

>30 


lO"*® 
16 - 23 


3x10® 
25 


10* 
150 


10® 10^ 
20 100 




2 


0.5 


1.1 


11.6 


(1) 


1x10^ - 3x10'' 
>40 












4 


0.4 


0.7 


8.5 


(1) 


1x10* - 6x10* 
20 - 60 




1x10® 
10 






N 166 


1 


0.3 


0.6 


10.5 


(2) 


5x10'^ - 2x10* 
25 - 150 












3 


0.3 


0.6 


12.6 


(2) 


3x10^ - 2x10* 
>30 












4 


0.4 


0.6 


16.2 


(2) 


4x10"^ - 2x10* 
>40 








10* 10^ - 10* 
20 30 - 60 


N 206 


1 


0.5 


0.8 


14.1 


(3) 


5x10^ - 3x10* 

>35 












2 


0.4 


0.5 


9.8 


(3) 


6x10^ - 2x10* 
20 - 80 










N 206D 


1 


0.3 


0.5 


4.8 


(3) 


1x10* - 3x10* 
10 - 20 










GMC 225 


1 


0.2 


0.4 


6.6 


(3) 


5x10^ - 2x10* 
15 - 20 












3 


0.3 


0.4 


6.7 


(3) 


7x10^ - 2x10* 
10 - 40 











Note. — Column (7); Results of LVG analysis Column (8)-(13); Results of previous studies, including both LVG and MEP analyses. 



References. — (1) Johansson et al. 1998; (2) Garay et al. 2002; (3) This work; (4) Heikkila ct al. 1999; (5) "Single-component fit" in Bolatto et al. 
2005; (6) "Cold Dense Component" and (7) "Hot Tenuous Component" of the "Dual-component fit" in Bolatto et al. 2005; (8) Kim 2006 
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Table6. Effect of X (CO) 







30 Dor No.l 


N 159 No.l 


N 159 No.4 


This work 










(Uniform X{CO)) 


X(CO) 


3x10-8 


3x10"*' 


3x10-*^ 




n(H2) (cm^'*) 


3x10-^ - 3x10^ 


3x10^ - 8x10^ 


IxlO^ - 6x10^ 




Tkin (K) 


>50 


>30 


20 - 60 


(Different X{CO)) 


X{CO) 


IxlO"** 


1x10"=^ 


3x10-'' 




n(H2) (cm-3) 


6xl03 - lxlO« 


2xl03 - 2x10^ 


lxl03 - 6xl03 




Tkin (K) 


>35 


>45 


20 - 60 


Heikkila et al. 1999 












X(CO) 


1.4x10-*^ 


1.2x10-5 


4.6xl0-'5 




n(H2) (cm-3) 


1x10^ 


3x10^5 


1x10^5 




Tkin (K) 


50 


25 


10 



-74- 



Table 7. Estimated PDR surface temperatures 



Region 


n 


Go 




References 




(cm-3) 




(K) 




30 Dor 


10* 


3500 


300 


1, 2, 3, 4 


N 159 


10* 


300 


100 


1,4,5 



Note. — Column(l): Regions. Column(2): Gas 
density. CoIumn{3): FUV flux in units of local inter- 
stellar value; 1.6x10^^ ergs cm~^ s~^. Column(4): 
Derived PDR surface temperature of the atomic gas. 
Column(5): References of FUV flux. 

^Estimation is done by using Figure 1 of Kaufman 
et al. 1999 

References. — (1) Bolatto ct al. 1999; (2) 
Poglitsch et al. 1995; (3) Werner et al. 1978; (4) 
Israel & Koornneef 1979; (5) Israel et al. 1996 



